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SUMMARY 

In-plane and interlaminar shear strength of a unidirectional SiC fiber-reinforced (BaSr)Al 2 Si 2 0 8 
celsian composite were measured by the double-notch shear test method between room temperature 
and 1200 °C. The interlaminar shear strength was lower than the in-plane shear strength at all tem- 
peratures. Stress analysis, using finite element modeling, indicated that shear stress concentration was 
not responsible for the observed difference in strength. Instead, the difference in layer architecture and 
thus, the favorable alignment of fiber-rich layers with the shear plane in the interlaminar specimens 
appears to be the reason for the low strength of this composite. A rapid decrease in strength was ob- 
served with temperature due to softening of the glassy phase in the material. 


I. INTRODUCTION 

Barium strontium aluminosilicate (BaSr)Al 2 Si 2 O g , commonly called celsian, is a potential matrix 
material for high temperature composites due to its low thermal expansion and excellent oxidation re- 
sistance (ref. 1). Its most useful crystalline form is monoclinic since it is phase stable up to about 
1600 °C (refs. 1 to 3). Like many other ceramics, the monolithic celsian is brittle and has relatively low 
flexural strength, -130 MPa. However, fiber reinforcement with e.g., 42 vol % of BN/SiC-coated con- 
tinuous Hi-Nicalon SiC fibers not only increases the flexural strength of the composite along the fiber 
direction drastically, -500 percent, but also, failure occurs in a non-catastrophic manner (ref. 3). 

Continuous fiber reinforced composites are orthotropic materials whose behavior lies between 
those of isotropic and anisotropic materials (ref. 4). Because the physical and mechanical properties 
are directional, the degree of this directionality strongly depends on the stacking sequence of fiber bun- 
dles. Therefore, the improvement to be realized in strength and failure behavior depends on how the 
composite is loaded with respect to principal material directions, in addition to the constituent properties 
and fiber/matrix interface bond strength. For example, since composites typically have low shear 
strength, shear failure could be dominant in the off-axis loading. Depending on composite architec- 
ture, such shear failures could take place either within the plane of composite, i.e, in-plane, or within 
the plane of thickness, i.e., interlaminar (ref. 5). These orientations are illustrated in figure 1. There- 
fore, to avoid unexpected failure and predict it beforehand by appropriate failure criteria, the factors 
affecting both shear strengths must be studied. 

A number of test methods have been developed to study both in-plane and interlaminar shear prop- 
erties. In-plane shear strength could be measured by an off-axis tension test (ref. 6), Iosipescu test 
(refs. 7 to 11), torsion test (ref. 12), etc. Interlaminar shear strength could be determined by short- 
beam test (refs. 13 and 14), Iosipescu test (refs. 10, 11, and 15), double-notch method (refs. 16 to 20), 
etc. The shear test method should lead to a uniform shear stress distribution on the shear plane. For 
example, a torsion test, which provides a "true" pure shear stress state in thin walled-tubes, could not 
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be used for plate composites. Even if it were used, processing ceramic tube test specimens can be pro- 
hibitively expensive. The Iosipescu test, which leads to a "nearly" pure-shear state, and was standard- 
ized by American Society for Testing and Materials (ASTM) for both polymeric (ref. 21) and ceramic 
composite materials (ref. 22), requires larger specimens. It is not suitable for elevated temperature 
testing with its present commercial metallic fixture. 

In the material development stage it is common to use test methods which are simple and require 
small specimens, although they may not provide the desired uniform stress-state. The double-notched 
shear (DNS) test, where notches make the shear stress distribution in the specimen highly nonuniform 
(refs. 16 to 20), was used in this study. Since the DNS test is inexpensive, the ASTM standardized it for 
both plastic (ref. 23) and ceramic matrix composites (ref. 22). The research community uses it fre- 
quently for comparative studies and quality checks. The DNS test was chosen in the present study be- 
cause it allows the measurement of both in-plane and interlaminar shear strength by the same test 
specimen configuration. Moreover, the DNS test could be used at elevated temperatures. 

Since the dimensions of the test specimens in this study were different from those in the ASTM 
Standard (ref. 22), and used by others (refs. 16 to 20), the first objective was to analyze the test speci- 
men by finite element modeling (FEM). In light of the findings from stress analysis, the second objec- 
tive was to carry out experiments up to 1200 °C to determine if there is a difference between in-plane 
and interlaminar shear strength of the celsian composite. 


II. STRESS ANALYSIS 

Composite materials are complex with several elastic constants due to their orthotropic structure. 
The stress analysis requires constants which may not be readily available, particularly when the speci- 
men material is new and made in small quantities. Fortunately, in such cases theoretical models could 
be used as alternative methods to estimate these constants from the constituent material properties. 
This approach, however, requires microstructural information including the number of constituents, 
constituent elastic constants, volume fractions, stacking sequence, and architecture of layers. 

Figures 1(a) and (b), respectively, show the interlaminar and in-plane shear specimens which are 
related by 90° rotation. The notches were made on the opposing faces perpendicular to the longitudinal 
axis, which is parallel to the fibers. The principal material directions were represented by the Cartesian 
coordinate system of 1, 2 and 3. As illustrated in figure 1, direction- 1 is along the fibers while direc- 
tions-2 and -3 are parallel to and perpendicular to the layers, respectively. Thus, extensional macro- 
scopic engineering constants E u , E 22 and E 33 represent the material coordinate system. Figures 2(a) 
and (b) show optical micrographs of the current composite perpendicular to and parallel to the fiber 
axis, respectively. It is clear that the composite microstructure contains distinct regions made up of pe- 
riodic "fiber-rich" and "matrix" layers. Although such a layered microstructure is unusual in unidi- 
rectional composites, it should be taken into account during stress analysis and evaluation of experi- 
mental strength values. 


II.l. Determination of Macroscopic Elastic Constants 

In the interlaminar specimen (fig. 1(a)), the layers are perpendicular to the plane of the figure, 
i.e., parallel to imaginary shear plane between the notches. In the in-plane specimen (fig. 2(b)), how- 
ever, the layers are in the plane of the figure, i.e., perpendicular to the shear plane between the notches. 
The first step in calculations was to determine the elastic constants of matrix and fiber-rich layers. The 
second step was to carry out further calculations based on findings in the first step and the number of 
layers present in the composite to obtain macroscopic elastic constants. 

The polycrystalline celsian matrix layer is isotropic and its two elastic constants are already known, 
so no additional calculations were necessary for it. For the unidirectional fiber-rich layer (laminae), 
however, in-plane elastic constants E L , E r , x > LT , and G LT were needed (L and T stand for longitudinal 
and transverse, respectively). For this, the constituent properties in table I were used. Since the fiber- 
rich layer had duplex coated fibers (nominal coatings thickness 0.4 pm BN and 0.1 pm SiC), in 
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addition to the matrix, the elastic constant calculations were made in three steps; first, for the BN 
coating on the Hi-Nicalon fiber, then for the SiC coating on the BN coated fiber and finally, for the 
celsian matrix surrounding the duplex coated fiber. The calculations were made such that the values 
obtained in each step were used as input for the following step, while taking into account the volume 
fractions of the constituents involved. For example, in the final step of the calculation of the longitudi- 
nal elastic modulus for the fiber-rich layer, E L , which was made for the celsian matrix, used the longi- 
tudinal modulus value obtained from the earlier step involving the SiC coating, E L/ , and was carried 
out using the rule-of-mixtures in equation (1). 


E l -E L fVf+E m (\-Vf) 


( 1 ) 


where E m is the elastic modulus of matrix and V f is the fiber volume fraction in the fiber-rich layer, 
which was about 66 percent. Similarly, the calculation of the transverse elastic modulus for the fiber- 
rich layer, E T , utilized the transverse fiber modulus obtained in the previous calculation step involving 
the SiC coating, £ r f , and was done using the Halpin-Tsai equation 1 (ref. 4): 


e t 


= E 


' 1 +&v/ 

l 1 -^/ J 


( 2 ) 


where 


(Er.flE m )- 1 

(E T j /£,„) + £ 


£ is equal to 2 for circular fibers. Similarly, the in-plane shear modulus of the fiber-rich layer, G LTy was 
calculated by equation (2) where E values were replaced by corresponding G values and £ is equal 
to 1. The in-plane Poisson ratio, v LT , was calculated by the rule-of-mixtures in equation (1). The re- 
sults of these calculations and volume fraction of the layers determined by the image analysis are 
shown in table II. 

For the interlaminar specimen (fig. 1(a)), the elastic constants needed for two-dimensional stress 
analysis were E u , E } 3 , G 13 and \) 13 . Because of microstructural considerations in figure 2, elastic con- 
stants E u and o 13 were calculated by equation (1) while £ 33 and G 13 were obtained by equation (2) using 
data in table II. The results are shown in table III. It should be mentioned that the composite elastic 
modulus along the fibers, E u , was also measured experimentally by strain gages and found to be 
165 GPa. As can be seen in table III, the corresponding calculated value is 172 GPa, i.e., a discrepancy 
of <5 percent. Considering the experimental uncertainties involving the measurements of volume frac- 
tions and elastic modulus values of the constituents, the agreement between the experimental and 
calculated values is very good. 

For the in-plane specimen (fig. 1(b)) the relevant elastic constants were E n , £ 22 , G 12 and \) 12 . Since 
the layers are parallel to the plane of the figure it is best to use laminate theory to determine the macro- 
scopic elastic constants. Laminate analysis requires stiffness of the matrix, the thickness and number 
of each layer, and the rotation angle between the principal material direction of each layer and the test 
specimen (refs. 4, and 24). The components of stiffness of the matrix for both layers were already 
determined in table II. Since the fibers are along the test specimen axis, the rotation angle for the fiber- 
rich layer is 0°. For the matrix, the rotation angle is not an issue because of isotropy. The thickness of 
the fiber-rich and matrix layers was determined to be 90±18 and 34±13 pm, respectively. Using this 


‘Halpin-Tsai equation is less accurate with fiber volume fraction approaching 1 (ref. 4), which was the case in calculations involving 
the coatings. However, since the calculations without coatings (only matrix and fibers) where fiber volume fraction was 0.42 also 
gave similar results, the degree of error may not be significant. 
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information, the so-called extensional stiffness matrix A(3x3) was determined. Macroscopic elastic 
constants for the in-plane (laminate) specimen were calculated by equation (3), (ref. 4). 


_ Al.l)A2,2) “^'(1.2) 
A(2,2)^ 

(3a) 

_ Al,l)A2,2) “^"(1,2) 
Au) A 

(3b) 

v _ A 1,2) 

12 - /i 

^2,2) 

(3c) 

^ A 6,6) 

^12 = . 

n 

(3d) 


where the indices in the parentheses show the location of components in the A matrix while h is the 
total specimen thickness. These results are also shown in table III. It should be mentioned that £,, cal- 
culated in equation (3a) agrees well with the value determined for the same constant by equation (1) in 
the interlaminar specimen. 

Based on figures 1 and 2 and constituent properties in table I, one would expect the elastic modulus 
to be highest along the fiber direction (£,,) and lowest perpendicular to the layers (£ 33 ), and intermedi- 
ate parallel to the layers (£ 22 ). In fact, the calculated values of E u (172 GPa), E 22 (153 GPa) and 
£33 (149 GPa) in table III are in agreement with this prediction. Notice, however, that the ratios of 
£ n /£ 22 , and £ n /£ 33 in this composite are not as large as those in the polymer matrix composites. 
Moreover, the layered microstructure in figure 2 did not much affect the ratio of transverse modulus 
values (£ 22 /£ 33 = 1 in the transversely isotropic composites). These results are believed to be due to the 
fact that the matrix elastic constants in the present composite are much larger and thus, a balance is 
present between the elastic constants and volume fractions of the Hi-Nicalon fiber, BN and SiC coat- 
ings, and celsian matrix. However, the above calculations showed that the modulus of the compliant 
BN coating, which was added to improve the fracture toughness by crack deflection, strongly affect the 
transverse elastic (£ 22 and £ 33 ) and shear (G 12 , and G 13 ) moduli. Therefore, the possible cracks in this 
coating as well as debonding at the fiber-matrix interface, which were not considered in these calcula- 
tions, could significantly reduce these values, resulting in larger ratios of the longitudinal to transverse 
moduli. 


II.2. Modeling 

Linear-elastic FEM was carried out to determine the effect of specimen orientation on the stress- 
distribution between the notches. Two-dimensional plane-stress linear elastic analysis was carried out 
by a commercial package 2 using eight-node quadrilateral isoparametric elements with isotropic and 
orthotropic elastic properties. As boundary conditions, a uniform compressive load was applied along 
the specimen axis, and the specimen was allowed to deform laterally. The models utilized the real 
specimen dimensions. To capture actual notch geometry, the notches were modeled as semicircles. 
Much finer mesh density was used at high stress regions surrounding the notches to improve the 
accuracy of calculations at these locations. Figure 3(a) shows the overall finite element model while 
figures 3(b) and (c) illustrate comparison of the notch region in the actual specimen (interlaminar) 
and in the model, respectively. The elastic constants used in the interlaminar and in-plane specimen 


2 ANSYS 5.5 University version, ANSYS Inc., Canonsburg, PA 15317. 
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models were E n , E 33 , G 13 and o I3 , and £ n , E 22 , G 12 and v ]2 , respectively. Since specimen dimensions 
were the same, both models were loaded by the same uniform compressive stress along the fibers. 


H.3. Comparisons of Stress Distributions 

Figure 4 shows the comparison of normalized shear stress distribution between the notches for inter- 
laminar and in-plane specimens. The normalization was made by the average stress, x mr , between the 
notches, equation (4). 


X 


ave 


db 


(4) 


where P a is the applied force while b and d are the distance between notches and depth of the specimen, 
respectively. It is clear in figure 4 that shear stress distributions in both specimens are not uniform. 
There is a high shear stress concentration near the notches, which is about 5.6 and 5.5 for the in- 
plane and interlaminar specimens, respectively. At exact notch roots, however, the shear stress is zero, 
because at these points the stress changes sign due to asymmetry in the notch geometry. Between the 
notches the shear stress decreases rapidly. Very similar shear stress distributions in the interlaminar 
and in-plane specimens imply that any difference between the experimental strength values between 
these specimens should not be attributed to the difference between shear stress distributions, but to 
other factors. This point will be discussed in detail later. 

Figure 5 shows the distribution of stress components between the notches in the interlaminar speci- 
men’. The stress values were normalized by the applied uniform stress, a a . As can be seen, the stress 
concentration, c v , near the notch roots is close to five times the shear stress concentration, x xy . Stress 
concentration along the x-direction, g v , is similar to that of x Ay * 4 . Notice, however, that both o v and 
o y stress concentrations are compressive. Therefore, the macroscopic failure of these specimens is 
expected to occur in shear between the notches. 


III. EXPERIMENTAL PROCEDURE 

Hi-Nicalon*' fibers, with an average diameter of 14 pm, were used as reinforcement. Prior to proc- 
essing, the fiber surfaces were coated with ~0.4 pm BN followed by -0. 1 pm SiC by chemical vapor 
deposition. The precursor to the celsian matrix of 0.75BaO-0.25SrO-Al 2 O 3 -2SiO 2 was made by solid- 
state reaction. The precursor powder consisted of mainly SiO : and BaAl 2 0 4 with small amounts of 
Ba 2 Si0 4 , a-Al 2 0 3 and Ba 2 Sr 2 Al 2 0 7 . This powder was made into a slurry with an organic solvent with 
various additives (ref. 25). Tows of BN/SiC-coated fibers were impregnated with the matrix precursor 
by passing them through the slurry. The resulting prepreg tape was dried and cut into pieces. Unidi- 
rectional composites were prepared by tape lay up (20 plies) and warm pressing. Finally, dense com- 
posites were obtained by hot pressing under vacuum in a graphite die. X-ray diffraction showed that 
the precursor was fully converted into the desired monoclinic celsian phase through solid-state 
reaction. The fully dense composite was about 4 mm thick and had about 42 vd % fiber. 

Shear strength was measured using 4x4x20 mm DNS specimens (fig. 1). The specimen cross- 
section was purposely made square to make one-to-one comparison between the in-plane and inter- 
laminar shear strength values at a given test condition without the interference of possible specimen 
volume/shear area effect. Notches, which had 0.3 mm width, were extended exactly to the middle of the 
specimens, i.e., d/2 = 2 mm, to ensure that failure takes place on the plane between the notch tips. 
Notches were in equal distance from the ends, and the distance between them was about 5 mm. 


Stress distribution was very similar in the in-plane specimen. 

4 Actual sign of shear stress between the notches is negative, as shown in figure 5. For convenience, shear stress values were shown 
as positive in figure 4. 

‘Nippon Carbon Co., Japan. 
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Dimensions of DNS specimens in this study were different from those recommended by the ASTM 
standard (ref. 22) because the test specimen in the standard has a rectangular cross-section and would 
not be suitable to study both in-plane and interlaminar shear strength. Moreover, a support fixture 
recommended in the standard to prevent possible buckling of specimens under compressive force was 
not utilized since lateral clamping of specimens could lead to artificially high shear strength values 
(refs. 17 and 22). However, based on the recommendation of a previous study (ref. 16), the specimens 
were machined to have a height to thickness (h/d) ratio of 1.87 to minimize buckling (fig. 1). Buckling 
of specimens was checked by two strain gages which were mounted on the opposite lateral faces. 6 The 
difference in strain values from these two strain-gages was <3 percent, indicating that the buckling ef- 
fect on strength values would be small. Force in both specimens was applied along the longitudinal 
axis by an MTS machine in displacement control with a cross-head speed of 2 pm/s. Nine interlami- 
nar and eleven in-plane shear specimens were tested up to 1200 °C in flowing N 2 . Microstructural 
analysis was carried out by X-ray and transmission electron microscope (TEM). TEM specimens were 
prepared by conventional methods. An ion-beam thinner was used to obtain electron transparency in 
the specimens. Fractography was studied by scanning electron microscopy (SEM) and optical micros- 
copy. 


IV. RESULTS 
IV. 1. Microstructure 

Figure 6 shows a typical TEM micrograph of a fiber with BN and SiC duplex coating and celsian 
matrix. The BN coating, which is amorphous, contains four successive layers and has an average total 
thickness of 0.4 pm. The SiC coating, however, is crystalline and has an average thickness of 0.1 pm. 
The BN layer w r as applied on the fiber surface to produce a compliant layer on the fiber and promote 
crack deflection away from the fibers while SiC was used to (i) form a diffusion barrier between BN 
and matrix, and (ii) to protect the fiber from damage during handling and composite processing 
(ref. 3). The compilation of various data showed that the thickness of the coating along with the bond 
strength of the coating/matrix and coating/fiber interfaces plays a critical role in determining the inter- 
laminar shear strength of composites (ref. 26). The X-ray diffraction (XRD) pattern taken from the 
polished section of composite, which is not shown here, showed monoclinic celsian as the only crystal- 
line matrix phase, indicating that the thermodynamically stable monoclinic phase formed in situ from 
the mixed oxide precursor during hot pressing. The unstable and undesired hexacelsian phase was not 
present. 


IV. 2. Mechanical Properties 

Figure 7 shows the comparison of typical load versus displacement curves for the interlaminar and 
in-plane specimens at 1000 °C. There are three important points which need to be emphasized here. 
First, the in-plane shear strength was higher than the interlaminar shear strength. Second, the failure 
behavior was different: in the in-plane specimens the failure always took place in a catastrophic man- 
ner while in the interlaminar specimens it usually occurred in stages. Third, although the shear area 
between notches was the same in both specimens due to their square cross-sections, shear stiffness of 
the in-plane specimens was usually higher than that of the interlaminar specimens. These observa- 
tions were independent of the test temperature used, and buckling of specimens was not observed. 

The results of the shear strength tests of both specimens as a function of temperature are shown in 
figure 8. The shear strength values reported here are the average stress between notches, T mr , which 
were calculated at maximum force by equation (4). Although equation (4) does not provide the intrinsic 
strength of the material, due to non-uniform stress distribution between the notches, it was useful here 


Measurement Group, P-3500 unit was used. 
MTS System Co. Eden Prairie, MN, USA. 
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since the goal was to compare the strength values of these two specimens under identical testing condi- 
tions. It is clear in figure 8 that the in-plane shear strength is higher than the interlaminar shear 
strength regardless of temperature. Another important finding is that the shear strength decreased 
with temperature in both specimens, which appears to be pronounced above 600 °C. 

Fractography showed some important differences. Figures 9(a) and (b) are the typical SEM 
micrographs of interlaminar and in-plane shear specimens, respectively, tested at 25 °C. The examina- 
tion of notch regions in the interlaminar specimens showed that the failure always occurred within the 
fiber-rich region, a strong indication that the fiber-rich layer was weaker. Notice also in figure 9(a) that 
there is a region where the fibers bridge the shearing surfaces. This is because the fiber-rich layers 
were not always parallel to the specimen axis, or the shear plane between the notches, i.e., some 
misalignment was present. As a result, cracks from the notches could not link up on the same shear 
plane although the notch roots were exactly in the middle of the specimens. Since the fibers are 
stronger than the matrix, the bridging fibers prevented catastrophic failure and continued carrying 
some of the load beyond its peak value. In fact, this was the reason why the load/displacement curves in 
the interlaminar specimens exhibited steps (fig. 7). Such failure features were also reported in the in- 
terlaminar specimens of carbon/epoxy composites (ref. 27) 

Fracture features in the in-plane specimens were markedly different as the complete separation of 
the specimen across the shear plane always occurred. This is because the shear plane in these speci- 
mens is perpendicular to fiber-rich layers: thus, any misalignment between the fiber-rich regions and 
specimen axis was not critical. As a result, the failure took place catastrophically at maximum force 
(fig. 7). Moreover, the in-plane specimens contained secondary cracks around the primary crack. 
Secondary cracks are believed to have occurred by the sudden release of elastic strain energy at peak 
force. Since no such cracks were observed in the interlaminar specimens, this indicates that the strain 
energy stored in the in-plane specimens during loading was higher. The comparison of stored elastic 
energy in figure 7, which is proportional to the area under the loading curve, is in agreement with this 
assessment. Figure 10 shows an important failure mode in the interlaminar specimen and has a 
significant implication as to whether or not the stress concentration at the notches or the location of a 
critical flaw determines where the fracture takes place. The crack on the left originated at the tip of the 
notch at the bottom; however, the crack on the right did not start at the tip of the top notch but 
surprisingly away from it, perhaps at a weaker site associated with the fiber-rich region. Results from 
such specimens were not included in figure 8, although the values were not much different. 


IV.3. Effect of Test Temperature 

To determine why the shear strength decreased with temperature (fig. 8), the shear planes and mi- 
crostructure were studied. The comparison of fracture planes in the interlaminar shear specimens, 
which were tested at 25 and 1200 °C, are shown in figures 11(a) and (b), respectively. The failure in 
both specimens occurred in the fiber rich layers. Notice that the surface of the room temperature 
specimen (fig. 1 1(a)) is relatively clean and cracks between the fibers and matrix fragments are clearly 
visible. The surface of the high temperature specimen (fig. 11(b)), however, has a distinct phase, which 
wetted the shear plane and glued the matrix fragments to the fibers. 

Because shear strength also decreased in the in-plane specimens, their fracture surfaces were also 
studied. However, these observations were not conclusive because of the alternating layers, and thus, 
thin foils were prepared for TEM. A typical dark field TEM micrograph is shown in figure 12. As 
marked, the pockets of amorphous phase were found between the celsian grains. Similarly, such a 
phase was also observed near the fibers in the fiber-rich layers. Since the amorphous glassy phase 
typically softens at intermediate temperatures, failure takes place at a lower external force through 
sliding of matrix grains and fibers. Therefore, reduction in shear strength with temperature in figure 8 
is believed to have occurred due to the softening of this amorphous phase. Such strength degradation 
by amorphous phase in ceramic systems is common and was observed previously by others as well 
(ref. 28). The deformation of grains through dislocation plasticity is unexpected since the temperature 
range in this study was not high enough for the activation of dislocation motion. Depending on the 
differences between the thermal expansions of constituents, the residual stresses could also affect the 
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strength below the processing temperature. However, the residual stresses could not be responsible for 
the observed reduction in strength with temperature here since their magnitudes are reduced with 
temperature. The difference in strength reduction with temperature in figure 8 may be related to the 
amorphous phase content on the corresponding fracture planes of these two specimens. 

It should also be mentioned that although strength decreased substantially with temperature as a 
result of softening of the grain boundary phase, the load/displacement curves did not reveal substantial 
nonlinear deformation features. One possible reason could be that the strain rate used in this study is 
high for such an observation. 


V. DISCUSSION 

V.l. Effect of Anisotropy on Stress Distribution in DNS Specimens. 

Distribution of stresses in the DNS specimens was shown to be dependent on such variables as the 
notch separation (refs. 16 to 20 ), notch depth (ref. 20 ), specimen width (ref. 16), notch distance from 
the specimen ends (refs. 16 and 17) and the presence/absence of support fixture (refs. 17 and 22). 
Thus, it was recommended that experimental strength measurements involving nonstandard DNS 
specimens should be complemented by stress analysis (ref. 17). Stress analysis of current specimens 
showed that shear stress distributions are nonuniform and there are large stress concentrations at the 
notches (fig. 4). 

Another variable which affects stress distribution in the DNS specimen is anisotropy, although this 
was not very obvious in figure 4 due to a small difference between the elastic constants of in-plane and 
interlaminar specimens. To demonstrate the effect of anisotropy on the shear stress distribution and its 
possible implications, two transversely isotropic unidirectional composites were compared with an iso- 
tropic material. The degree of anisotropy in these composites is represented by E u /E 33 (or E u /E 22 ) 
(ref. 16). Graphite fiber reinforced (£,,/ £33 = 12.8) and glass fiber reinforced (E n /E 33 = 3.3) polymer 
composites were chosen because their elastic constants are well characterized and anisotropy between 
them is large enough for discussion here (ref. 29). Note that anisotropy in the current specimens was 
small in spite of its layered microstructure, i.e., EJE^ = 1.12 for the in-plane and E U IE 33 = 1.15 for 
the interlaminar specimens, and it was nearly transversely isotropic, i.e., E^JE^ ~ 1. For isotropic 
material, celsian was selected. The elastic constants used in the analysis were 
E u = 155 GPa, E 22 = E 33 = 12.1 GPa, G 12 = G l3 = 4.4 GPa and v n = u 13 = 0.248 for the graphite 
composite, and £,, = 50 GPa, E 22 = E 33 = 15.2 GPa, G , 2 = G 13 = 4.7 GPa and x> ]2 = \) 13 = 0.254 for the 
glass composite. Elastic constants for celsian were E = 96 GPa and x> = 0.32 (table III). 

Figure 13 compares the behavior of different types of materials. It is clear that anisotropy has a 
strong effect on the stress distribution, particularly on the stress concentration near the notches. Stress 
concentration was largest in celsian (~5.8), intermediate in the less anisotropic glass composite (~3.6) 
and smallest in the most anisotropic graphite composite (~2.3). Such findings and trends are in 
agreement with a similar stress analysis involving the S-Glass/Epoxy, Graphite/Epoxy and Graph- 
ite/Polyimide composites (ref. 16), where the E n /E 33 ratios of the composites were 2.4, 7.3 and 43.2, 
respectively. Notice in figure 4 that shear stress concentration in the present composite was ~5.6, which 
is not much different from that of celsian in figure 13. This analysis clearly shows that anisotropy of 
the unidirectional composite reduces the stress concentration at the notches, and as a result, makes the 
stress distribution much more uniform. The decrease in stress concentration with anisotropy in figure 
13 should not be surprising since the smaller lateral stiffness E 33 (or E 22 ) material 
increases the lateral deformation. The implication of figure 13 is that the DNS specimen produces 
better quality strength data if the test involves the composite with large anisotropy, assuming that the 
specimen is loaded along the stiffest£ n direction. The strong effect of anisotropy on the stress distribu- 
tion of composites was also shown in the losipescu specimen (ref. 10 ). 
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V.2. Interlaminar versus In-Plane Shear Strength 


Since the notch magnifies the effect of far-field applied force, the stress concentration associated 
with it could have a big impact on the measured strength values. In fact, such impact in DNS speci- 
mens was demonstrated in ceramic composites where the strength decreased with the magnitude of 
stress concentration as well (refs. 16 and 19). 

In this study, the stress concentration in both specimens was almost the same in spite of the layered 
microstructure. 


^interlaminar — ^f,in-plane (^) 

Based on stress concentration alone, one would have expected the interlaminar and in-plane shear 
strength to be similar. However, results in figure 8 showed that the interlaminar strength was lower 
than in-plane strength at all temperatures, i.e., 

^13 (interlaminar) < ^12 (in-plane) (6) 

This discrepancy between the prediction and the experiments indicates that the stress concentration 
could not be responsible for the difference observed in strength values. Moreover, even if the calculated 
elastic constants in table III had underestimated the "actual" anisotropy, 8 this conclusion would have 
still been valid since the anisotropy, which is higher in the interlaminar orientation, is expected to 
reduce the stress concentration and thus increase the strength of interlaminar specimens. Since this 
was not the case in the experiments, this clearly excludes the stress concentration as the possible reason 
for the difference in strength. 

Instead, other factors such as microstructural features including the size and location of the 
strength limiting cracks and pores, which are not included in the continuum stress analysis, are 
believed to be responsible for the strength difference. In the 0/90° composites, for example, the differ- 
ence between the in-plane and interlaminar shear strength was attributed to the presence/absence of 
interactions between the matrix crack with fibers (ref. 5). In the interlaminar specimen this interaction 
is weak while in the in-plane specimens it is strong; as a result, the in-plane shear strength is always 
higher than the interlaminar shear strength in 0/90° composites (ref. 5). The same types of interactions 
with fibers are not expected here because the present composite is unidirectional (no 90° fiber-bundles) 
and shear planes in both specimens are parallel to the fibers. 

As indicated earlier, the failure in the interlaminar specimens always occurred within the fiber-rich 
layer which is the multiphase region with a high density of defects and grain boundaries. So, the proc- 
ess induced defects such as cracks and pores are often found at the fiber/matrix interface where initia- 
tion and growth of delamination cracks could take place when the crack path is parallel to applied 
stress, such as shear forces (ref. 17). An important observation on the critical role of the defective fiber- 
rich region on the shear failure of an interlaminar specimen was already shown in figure 10. The fact 
that the crack initiated away from the notch in figure 10 clearly indicates that the location of the failure 
and the apparent strength of the composite is determined not by the high stress concentration near the 
notch, but by local stresses near the most critical defect in the microstructure. Because such defects are 
contained within the fiber-rich layer, strength is determined by the orientation of this layer with respect 
to the shear plane: when it is parallel and perpendicular, it leads to low (interlaminar) and high 
(in-plane) shear strength values, respectively. 

High defect concentration in the fiber-rich layers may also be indicated indirectly by the slope of the 
load/displacement curves in figure 7 where the interlaminar specimen is less stiff. Since stiffness is 
defined as displacement per unit load, this suggests that the ligament surrounding the shear plane in 
the interlaminar specimens contains a higher fraction of the compliant part in the composite, i.e., fiber- 
rich layers. Although the difference is small, the calculated shear modulus values in Table III also 


8 Equations used to calculate unknown elastic constants assume perfect bonding at the interfaces and no porosity /cracks at the interface and in 
the matrix. Therefore, actual constants are expected to be different from those in table III based on the deviation from the assumed case. 
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show a similar tendency, i.e., 58 GPa and 56 GPa for the in-plane and interlaminar specimens, 
respectively. 

Figure 5 showed that normal stress components c v and c y on the shear plane are compressive. 
However, in a real microstructure containing cracks and pores, both shear and tensile stress compo- 
nents could be present, 9 i.e., mixed mode. Assuming the presence of a small precrack within the fiber 
rich layer, the stress intensity values of interest are K u and AT,, i.e., Mode-II and Mode-I, respectively. 
With the condition that the crack does not deviate from the original plane, which is expected here since 
the crack plane is parallel to the fibers, the fracture condition for mixed-mode loading is given by 
(ref. 30): 



V K ttc ) 


Bg^ko^ 


K 


Ic 


= 1 


(7) 


where A and B are the geometric factors, a c is the critical crack size and K /c and K IIC are the critical val- 
ues of stress intensity. So, the failure is expected to initiate at a critical crack whenever the combination 
of shear stress x and tensile stress a satisfy the above condition. Note that regardless of mode-mixity, 
equation (6) shows that stress to failure (strength) is inversely proportional to crack length, which 
means that the material with the largest initial crack would likely have lower strength. Therefore, the 
difference between shear strength of in-plane and interlaminar specimens could be attributed to the 
difference in the defect density on the corresponding shear planes in these two specimens, i.e., lower 
interlaminar shear strength is due to favorable orientation of the defective fiber-rich region with the 
shear plane. The area/volume could also have a big impact on the measured strength values of brittle 
materials. However, this is not expected to be important in this study since shear areas in both speci- 
mens were the same. 


VI. SUMMARY AND CONCLUSIONS 

Stress analysis showed that the shear stress distributions in interlaminar and in-plane specimens 
were non-uniform and a large stress concentration was present at the notches. However, stress distri- 
butions in both specimens were almost the same due to similarities in the corresponding elastic 
constants. It was also shown that anisotropy of the unidirectional composite is important in stress 
distribution and reduces the stress concentration at the notches. The reason for high stress concentra- 
tion in current specimens was due to relatively high E 22 and E n . 

Interlaminar shear strength was lower than in-plane shear strength at all temperatures. Moreover, 
the interlaminar specimens failed in stages while die in-plane specimens fractured catastrophically at 
maximum load. This could be attributed to the difference in the layer architecture. 

Comparison of results from stress analysis and experiments suggests that the difference between 
the strength in the two orientations cannot be explained by the differences in stress concentration, but 
rather is due to the favorable alignment of highly defective fiber- rich layers with the shear plane in the 
interlaminar specimens. Since these layers were perpendicular to shear plane, the strength of in-plane 
specimens was higher. 

The strength of both in-plane and interlaminar specimens decreased with temperature above 
600 °C. Microstructure analyses suggest that this was due to softening of the amorphous glassy phase 
in the material. 


9 Any failure away from the plane between notches, such as the one in figure 10, have more tensile component. 
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TABLE I.— CONSTITUENT ELASTIC PROPERTIES 
AND VOLUME FRACTIONS FOR THE HI-NICALON 
FIBER-REINFORCED CELSIAN COMPOSITE 


Constituent 

Elastic 

modulus, 

GPa 

Poisson's 

ratio 

Volume fraction, 
percent 

Hi-Nicalon fiber 

270 

0.20 

42 

Celsian matrix 

96 

0.32 

51.8 

BN 

62 

0.17 

4.9 

SiC 

450 

0.21 

1.3 


TABLE H.— CALCULATED IN-PLANE 
ELASTIC CONSTANTS FOR THE LAYERS 
IN FIGURE 2. THE SUBSCRIPTS L AND T 


INDICATE THE LONGITUDINAL AND 
TRANSVERSE DIRECTIONS IN THE 
MATERIAL COORDINATE SYSTEM 


Property 

Fiber-rich 

layer 

Matrix layer 

E, (GPa) 

200 

% 

Er (GPa) 

174 

96 

^LT. 

0.24, 0.21 

0.32 

G IT (GPa) 

66 

36 

Volume fraction 

0.73 

0.27 


TABLE III. — CALCULATED 
ELASTIC CONSTANTS FOR 
THE ENTIRE COMPOSITE 


Elastic constant 

Value 

E„ 

172 GPa 

E„ 

153 GPa 

E» 

149 GPa 


56 GPa 

G i: 

58 GPa 


0.26 

_ai 

0.25 
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Figure 1 . — Double notched shear (DNS) test specimens used in the 
current study: (a) interlaminar and (b) in-plane. Notice that they are 
related by 90°. 
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Figure 3. — Finite element model (a) of the DNS test specimen, and comparison 
of the notch geometry in the actual interlaminar specimen (b) and model (c). 


N AS A/TM — 2000-2 1 0608 


15 


Normalized stress, a x /o- a , or y /a a , T xy /cr a Normalized shear stress, T xy /T aver age 


6 



Normalized distance between notches 


Figure 4. — Comparison of normalized shear stress distribution between the 
notches for interlaminar and in-plane specimens. 



Figure 5. — Distribution of relevant stress components between notches in the 
interlaminar specimen. 
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Figure 6. — TEM micrograph of interface region between Hi-Nicalon fiber and celsian 
matrix showing the duplex BN and SiC coatings. Notice that the BN coating con- 
sists of four layers. 



Celsian 


Hi-Nicalon 

fiber 



Figure 7. — Comparison of typical load versus displacement curves of the inter- 
laminar and in-plane specimens at 1000 °C. 
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Figure 8. — Temperature dependence of the shear strength of interlaminar and in- 
plane specimens. 
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Figure 9. — SEM micrographs showing typical shear fracture features in the (a) interlaminar and (b) in- plane 
specimens tested at room temperature. 



Figure 10. — SEM micrograph showing a failed 
interlaminar specimen. Initiation of an unusual 
shear fracture away from the notch indicates that 
the location of failure is determined by the local 
stress-state associated with the critical defect in 
the fiber-rich regions rather than the stress concen- 
tration near the notches. 
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Figure 1 1 . — SEM micrographs showing fracture 
surfaces of interlaminar specimens tested at 
(a) room temperature and (b) 1 200 °C. Notice 
the absence and presence of a viscous phase 
in (a) and (b), respectively. 



Figure 12. — Dark field TEM image showing the 
amorphous grain boundary phase between 
the grains of celsian matrix. 
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Figure 13. — Effect of material anisotropy on the shear stress distribution between 
the notches. 
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